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Abstract 

We discuss several experimental and theoretical techniques historically used for 
Hasselmann equation wind input terms derivation. We show that recently de¬ 
veloped ZRP technique in conjunction with high-frequency damping without 
spectral peak dissipation allows to reproduce more than a dozen of fetch-limited 
field experiments. Numerical simulation of the same Cauchy problem for differ¬ 
ent wind input terms has been performed to discuss nonlinearity implications 
as well as correspondence to theoretical predictions. 
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1. Introduction 


It is generally accepted nowadays that ocean surface wave turbulence is de¬ 
scribed by Hasselmann equation (hereafter HE) 


de 

dt 


duly de 


= S, 


nl 


Sin 


Jdis 


dk dr 

for spectral energy density e = e(k,f,t), wave dispersion to = 
ear, wind input and wave-breaking dissipation terms S n i , Si n 


(1) 

6i{k) and nonlin- 
and Sdiss corre- 
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spondingly. 

While this acceptance implicitly assumes that HE is some sort of mathemat¬ 
ical reduction of primordial Euler equations for incompressible fluid with free 
surface, it is formally true, in fact, only for advection and four-wave 


interaction S n i terms. 

As far as concerns St n and Sdiss terms, there is no consensus in the worldwide 
oceanographic community about their parameterization. To our belief, it is one 
of the reasons, indeed, for “tuning knobs” (adjusting coefficients) necessity in 
operational models for their adjustment to different ocean situations set-ups. 
Another reason for using “tuning knobs” is the underestimation of the leading 
role of S n i■ In other words, the role of the “tuning knobs” also consists in 
“undoing” the deformation incurred to the primordial equations model through 
substitution of the exact nonlinear term S n i with DIA -like simplifications. 

We believe that such currently widely accepted approach of using ” tuning knobs” 
is conceptually misleading, and continuing efforts on improving the source terms 


fl, fl¬ 


are fruitless because the nonlinear term S n i is the leading term of HE 
All other source terms are, in a sense, relatively small corrections in significant 
ranges of frequency. 

Dominations of the nonlinearity exhibit itself in HE in the form of self-similar 
solutions, observed in the field and numerical experiments fl, fl fl fl fl fl- 
Self-similarity analysis in conjunction with the field and numerical experiments 
analysis allowed to build new ZRP wind input term fl - analytical solution of 


HE. 


In the current paper we use alternative approach to HE simulation, which in 
addition to new ZRP wind input term uses another physically based assump¬ 
tions - absence of the spectral peak dissipation and implicit wave damping due 
to wave-breaking. We justify this new approach through comparison with more 
than the dozen of field measurements. 

Another result of current paper is the development of the set of tests, based 
on self-similarity analysis, allowing to make the judgment relative to correct¬ 
ness of arbitrary numerical simulation performed in the frame of HE through 
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comparison of the observed characteristics of wave ensemble with theoretical 
prediction. 


2. Current state of wind input source terms 


Nowadays, the number of existing models of Si n is large, but neither of them 
have firm theoretical justification. Different theoretical approaches argue with 


each other. Detailed description of this discussion can be found in the mono¬ 
graphs and the papers 0 , 0, 0, 0, 0. 


The development of wind excited waves models has begun as far back as 20-ies 
of the last century in the well-known works of Jeffreys 0|, 0|. His model is 
semi-empiric and includes unknown ” sheltering coefficient”. All other existing 
theoretical models are also semi-empiric, with one exclusion famous Miles 
model 0. This model is rigorous, but is related to idealized situation - initial 
stage of waves excitation by laminar wind with specific wind profile U(z). 

Miles theory application is hampered by two circumstances. First is the fact 
that atmospheric boundary layer is the turbulent one, and creation of rigorous 
analytical theory of such turbulence is nowadays unsolvable problem. 

There is the opinion, however, that wind speed turbulent pulsations are small 
with respect to horizontal velocity U(z), and they can be neglected. This doesn’t 
mean that turbulence is not taken into account at all. It is suggested that the 
role of the turbulence consists in formation of the averaged horizontal velocity 
profile. 

This widely spread opinion is that horizontal velocity profile is distributed by 
logarithmic law 

U(z) = 2.5u* In— (2) 

zo 

Here it* is friction velocity and zq - the roughness parameter 

,2 


uZ 


zo = C ch S- 
9 


( 3 ) 


where C c h ~ 10 2 is experimental dimensionless Charnok constant. 
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One should note that appearance of anomalously small constants, not having 
’’formal justification”, is extremely rare phenomenon in physics. Eq. ©, © 
mean that roughness parameter is very small: for typical ocean conditions - 
wind speed 10 m/sec on the height z = 10 m we get z o — 10~ 4 m. Such 
roughness is only twice the size of viscid layer, defined from multiple experiments 
on turbulent wind flow over smooth metal plates. 

Usage of Eqs. ©, © assumes therefore that ocean behaves as smooth metal 
surface. This is not correct. Horizontal momentum is transferred to the smooth 
plate on it’s surface itself, while in the ocean this process happens differently. 
Momentum offtake from atmospheric boundary layer is smoothly distributed 
over the whole width of the boundary layer and begins from the highest ’’con¬ 
currence layer”, i.e. from the height where phase speed of the fastest wave 
matches the horizontal velocity. 

Momentum offtake leads to horizontal velocity distribution U{z) dependence on 
time, waves development level and energy spectrum. Meanwhile, Miles insta¬ 
bility increment is extremely sensitive to the horizontal velocity profile (there is 
no waves excitation for linear profile U(z) for Miles theory, for example). The 
velocity profile is especially important for slight elevations of the order of several 
centimeters over the water surface, which is almost unknown and difficult for 
experimental measurements. However, there are some advances in this direction 
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velocity 


The necessity of taking into account of waves feedback on the horizontal 
profile has been understood long time ago in the works of Fabrikant 
Nikolaeva, Zyrling 22] later continued in the works of Jannsen H and explained 


in details in the monograph [llj in the form of ”quasi-laminar” theory. This 
theory is not accomplished. 

To consider the theory as self-consistent even in the approximation of turbulence 
absence, it is necessary to solve equations describing horizontal velocity profile 
U ( 2 ) together with Hasselmann equation, describing energy spectrum evolution. 
This is not done yet. 

Aside that fact, many theoreticians do not share share the opinion about tur- 


211 and 


4 







bulent pulsations insignificance, and consider them as the leading factor. Cor¬ 


responding TBH theory by Townsend, Belcher and Hunt m is alternative to 
quasi-laminar theory. Both theories are discussed in H- 

There is another approach, not connected with experimental analysis - numeri¬ 
cal simulation of boundary atmospheric layer in the frame of empiric theories of 
turbulence. It was developed in the works |l3, Q, H, Si- Since those theories 
are insufficiently substantiated, the same relates to derived wind input terms. 
For all the variety of theoretical approaches of Si n definitions, all of them are 
are ” quasi-linear”, assuming: 


Sin = 7(w, </>)e(w, <j>) 

where standard relation 


7(w,</>) = —uP(—,4>) 

Pw Mp 

is being used. Here co p = -, where u is the wind speed, defined differently in 
individual models. Function /3 is dimensionless and is growing with the growth 
of -f. 

However, even for the most ’’aggressive” models of wind input terms the value 
of /3 does not exceed several units, but usually 0 < /? < 1. In some models 
(see,for example 


15|) /3 becomes negative for the waves propagating faster than 


the wind, or under large angle with respect to the wind. 

Looking at multiple attempts of Si n experimental definition, one should note 
that all of them should be carefully critically analyzed. That criticism is not 
about the integrity of measurements itself, but about the used methodology and 
data interpretation correctness, and the possibility of transfer of the conclusions 
made in artificially created environment to real ocean conditions. 

Another significant amount of experiments, belonging to so-called ’’fractional 
growth method” category, has been performed through energy spectrum mea¬ 
surement in time and calculation of the corresponding 7 through 
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(4) 


7(w,</>) 


1 de(ui, (/)) 
e(u>, (j>) dt 


Eq.® is, in fact, the linear part, or just two terms of the HE Eq.JTJ. This 
method is intrinsically wrong, since it assumes that either advection and 

nonlinear S n i terms of Eq. (ID are absent together, or relation 


dojde_ 

~d%d? = ' 


(5) 


is fulfilled. 

First assumption is simply not correct, since neglected terms are defining in 
ocean conditions. The second assumption is also wrong since Eq.® contradicts 
Eq.®. Therefore, in the relation to ’’fractional growth method” we are just 
citing the single relevant publication by Plant Q where, it seems, author well 
understood the scarcity of the ’’fractional growth method”. 

As a matter of fact, the real interest present the experiments, which used mea¬ 
surements of the correlation between the speed of the surface growth and the 
pressure to the surface: 


Q(uj) = Re < rj(uj)P*(u]) > 


( 6 ) 


Unfortunately, the number of such experiments is limited, and not all of them 
have significant value for ocean phenomena description. Also, one should take 
out of consideration the experiments performed in laboratory conditions. 
Consider, for example, the set of experiments described in 27]. These exper¬ 


iment were performed in the wave tank of 40 m length and 1 m depth. Ex- 
perimentators created the wind blowing at the speed up to 16 m/sec, but they 
studied only short waves no longer than 3 m, moving no faster than 1.3 m/sec. 
Therefore, they studied the very short-wave tail of the function /3 in the condi¬ 
tions far from the ocean ones. The value of these measurements is not significant. 
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The same arguments relate to multiple and precisely performed measurements in 
the Lake George, Australia 28]. The depth of this lake, in average, is about 1 m. 


That is why on its surface can propagate the waves not faster than 3.3 ml sec. 
The typical wind speed, corresponding these measurements, was 8 — 12 ml sec. 
Therefore, while the results of these measurements are quite interesting and 
correspond to theoretical predictions [29(, obtained expression for Si n is quite 


arguable, not only because is non-improvement to ’’quasi-linear” theory, but 
also being in complete contradiction with it. 

After critical analysis of experiments on Si n measurements, only three of those 


deserve an attention. Those are the experiments by Snyder et al 
Shemdin 


• [.30], 


Hsiao, 


311] and Hasselmann, Bosenberg [32]. The experiments were performed 


in the open ocean and measured direct correlations of surface speed change and 
the pressure. 

Those experiments were performed fairly long time ago, their accuracy is not 
quite high and scatter of data is significant. Therefore, their interpretation is 
quite ambiguous. Anyhow, these experiments produced two well-known formula 
for f3. For Snyder and Hasselmann, Bosenberg experiments: 


P = 0.24(£ — 1), £=—cos</> (7) 

w 0 

P = 0, £<1 (8) 

and for Hsiao-Shemdin 

P = 0.12(0.85£ - l) 2 (9) 

P = 0 otherwise (10) 

The difference between these S^ n formula is significant. Comparison of wind 
forcing performed on measured spectra Q shows that Snyder-Hasselmann-Bosenberg 
form gives 5 4- 6 times bigger value of S)n than Hsiao-Shemdin one. 

Furthermore, the Hsiao-Shemdin form agrees with Jeffreys theoretical model, 
while Snyder-Hasselmann-Bosenberg one is in disagreement with any known 


7 






theoretical models. 

Summing up, we can conclude that at the moment there is no solid parameter¬ 
ization of Si n accepted by worldwide oceanographic community. Keeping that 
fact in a mind, we decided to go our own way - not to build new theoretical 
models and not to reconsider both old and completely new measurements of 


Sin ■ 

For almost seventy year, counting from works of Garrett and Munk 


33|, physi¬ 


cal oceanography assimilated tremendous amount of experimental facts on basic 
wind-wave characteristics - wave energy and spectral peak frequency as a func¬ 
tions of limited fetch. Such experiments are described in works Q, 0, 7, Q]. 
From the other side, numerical methods of Hasselmann equation [1] solution for 
exact term S n i and specified in advance terms Si n and S<u S s have been improved 
significantly. This can be done not only for duration-limited domain, but for 
fetch-limited domain too. 

Therefore, we proposed purely new pragmatic approach to definition of Si n . We 
have chosen Si n function in such a way that numerical solution of Hasselmann 
equation explains maximum amount of known field experiments.The result was 


the Sin function described in details in [^land named thereafter ZRP function. 
It is important to emphasize that work 4] assumed localization of energy dis¬ 
sipation in short waves. This assumption contradicts widely accepted concept, 
but we explain the difference in the following chapter. 


3. Two scenarios of wave-breaking dissipation term: spectral peak or 
high-frequency domination? 

In current section we explain why there is no need to use dissipation in the 
spectral peak area. 

The spectral peak frequency damping is widely accepted practice, and is in¬ 
cluded as an option in the operational models WAM , SWAN and WW3. His¬ 
torically, it was done apparently by need. 

The necessity was caused by wind input function Si n in Snyder form. Fast wave 





energy growth was observed in no-dissipation calculations, which didn’t match 
results of field measurements. Despite the result was obtained with the help of 
DIA model of S n i, it is qualitatively correct, because is also confirmed by our 
numerical calculations using exact nonlinear interaction term S n i. 

It is shown below that Snyder wind input without dissipation gives 5 — 6 times 
bigger energy growth than other tested wind input functions ( ZRP , Tolman- 
Chalikov and Hsiao-Shemdin). This doesn’t mean, however, that long-wave 
dissipation exist, indeed. From our viewpoint, the necessity of its introduc¬ 
tion is explained by Snyder model imperfection, based on not quite accurate 
experiments. 

We don’t see the physical reasons for energy-containing long waves breakings. 
Their steepness in the conditions of typically developed wave turbulence is not 
big: g =< V?y 2 > 1|/2 ~ 0.1, or even smaller. Because this value is very far from 
limiting steepness of Stokes wave gs — 0.3, thise waves are essentially weakly- 
nonlinear. Besides those waves, more shorter waves inevitably develop, having 
the steepness approaching to the critical one, and those waves break. There is 
nevertheless no reason to expect that these waves have the same phase velocity 
as the energy-containing ones. 

Unfortunately, the theory of ’’wave-breakers” is not developed yet. In our view, 
which we don’t consider based enough, one of the possible variants of such theory 
could be the following. 

The primordial Euler equations for potential flow of deep fluid with free surface 
has the self-similar solution 


r,{x,t)=gt 2 F( — 


( 11 ) 


This solution was studied numerically in the framework of_simplified MMT 
(Maida — McLaughlin — Tabak) model of Euler equations 34]. 

In Fourier space this solution describes the propagation to high wave-numbers 
and returning back to dominant wave spectral peak of fat spectral energy tail, 
corresponding in real space to sharp wedge formation at time t = 0 and space 
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point x = 0. This solution describes formation of the ’’breaker”. 

In the absence of dissipation, this event is invertible in time. Presence of high- 
frequency dissipation chops off the end of the tail, just like “cigar cutter”, and 
violates the tail invertability. Low and high harmonics, however, are strongly 
coupled in this event due to strong nonlinear non-local interaction, and deformed 
high wave-numbers tail is almost immediately returns to the area of spectral 
peak. As soon as fat spectral tail return to the area of the spectral peak, total 
energy in the spectrum diminishes, which causes settling of the spectral peak 
at lower level of energy. This process of ” shooting” of the spectral tail toward 
high wave-numbers, and its returning back due to wave breaking is the real 
reason of ’’sagging down” of the energy profile in the spectral peak area, but 
was erroneously associated with the presence of the damping in the area of 
spectral peak. 

This explanation shows that individual wave-breakings studies 
the proof of spectral peak damping presence. 

Also there is another, direct proof of the fact that the damping is localized 
in the area of short waves. It is the measurements of quasi-one-dimensional 
’’breakers” speed propagation - strips of foam, which accompany any developed 
wave turbulence. Those airplane experiments, recently performed by P.Hwang 
and his team 


29|, [351 are not 
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38U39(, show that wave breakers propagate 4-5 times slower 


than crests of leading waves. 

Based on the above discussion, we propose to use only high-frequency damping 
as a basis of alternative framework of HE simulation. One can “implicitly” 
insert this damping very easily without knowing its analytic form via spectral 
tail continuation by Phillips law ~ w -5 . 

Replacement of high-frequency spectrum part by Phillips law is not our inven¬ 
tion. It is the standard tool offered as an option in operational wave forecasting 
models, known as the ’’parametric tail”, and corresponds to high-frequency dis¬ 
sipation, indeed. For the practical definition of Phillips tail it’s necessary to 
know two more parameters: coefficient in front of it and starting frequency. 
The coefficient in front of w~ 5 is not exactly known, but is unnecessary to be 
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defined in the explicit form - it is dynamically determined from the continuity 
condition of the spectrum. As far as concerns another unknown parameter - 


the frequency where Phillips spectrum starts - we define it as fo = ^ = 1.1 Hz 
as per Resio and Long experimental observations 


40 


111 ]. 


That is the way the high frequency implicit damping is incorporated into al¬ 
ternative computational framework of HE. We think that the question of finer 
details of high-frequency damping structure is of secondary importance at cur¬ 
rent stage of alternative framework development. 


4. Checking of the new modeling framework against theoretical pre¬ 
dictions and field measurements 

To check alternative framework for HE simulation, we performed numerical 
tests for waves excitation in limited fetch conditions. As it was already men¬ 
tioned, alternative framework is based on exact nonlinear term S n i in WRT 
form and ZRP new wind input term: 


Sin = 7 e 

( 12 ) 

/ \ 4 / 3 

7 - 0.05 Pair u W m 

(13) 

Pwater \^0 / 


f cos 2 6 for — 7t/2 < 9 < 7t/2 


m = { 

(14) 

1 0 otherwise 


9 Pair 1 O 1 n—3 

luq =- , -= 1.3 • 10 

(15) 


^10 Pwater 


The coefficient 0.05 in front of Ea. lflTil) was found through carefully performed 
numerical experiments with different coefficient values to get the best corre¬ 
spondence with experimental data. 

Thus 


j3 = 0.05w 


-) 

w 0 ) 


4/3 


m 


(16) 


We have chosen /3 in the form of power function of frequency for the following 
reasons. It is well known from various field experiments [ 6 i] that wave energy and 
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spectral frequency maximum dependencies on the fetch are the power functions 
of fetch: 


= £oX p 

(17) 

= woX -9 

(18) 


This observation is in excellent agreement with the fact that conservative sta¬ 
tionary Hasselmann equation 


dui de 

at r \ ^ Til 

Ok ox 

has two-parameter family of self-similar solutions 



6 , 



(19) 


e = X P+q F(u X q ) (20) 

which lead to dependencies Eqs. (Ha, m- 

One of the p , q is free, but they are connected by relation 


10p-2q = l (21) 

called thereafter the ” magic law”. 

Analysis of field experiments Q shows that ’’magic law” is fulfilled with high 
accuracy in many of them. 

HE Eq.J3 has self-similar solutions only if it is supplied by wind input term in 
power form: 

P = 0J s f(6) (22) 

and self-similar substitution gives in this case q = 2^7 • the mos t known 
experiments p = 1 and q = 3/10. That gives us the value s = 4/3. 

It is important to note that if is not the power function, the Eq.([3 

solution in the absence of long-wave dissipation is ” quazi-self-similar” in typical 
cases. In this case p and q are slow functions of the fetch, but the ’’magic 
law” Ea. (l2Tl) is still fulfilled Q. Strictly speaking, an existence of the universal 
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ZRP forcing , wind speed U=10 m/sec 




xg/U A 2 


(a) Solid line - numerical experiment, (b) Exponent p of the energy growth 
dashed line fit by 2.9 • 10 -7 • as the function of fetch x 


Figure 1 


for any ocean conditions expression of /3 {to, 9) is not proved, because turbulent 
boundary layer is different for typical oceans, passats (trade winds) or mountain 
coastal line. Therefore, the fact of explanation by ZRP expression for of at 
least half of known field experiments can be considered as big success. 

FigH] shows that total energy is growing along the fetch by power law in accor¬ 
dance with Eq. El with p = 1.0. 

Dependence of mean frequency on the fetch, shown on Figj2] also demonstrates 
perfect correspondence of numerical results and corresponding self-similar de¬ 
pendence Ea. (fl8l) with q = 0.3. 

Fig l3al presents directional spectrum as a function of frequency in logarithmic 
coordinates. One can see that energy curve on the left figure consists of segments 
of: 


1. Spectral maximum area 

2. Kolmogorov-Zakharov spectrum 

3. Phillips high frequency tail w -5 

Figl3bl shows log-log derivative of the spectral curve from Fig|3a| figure, which 
corresponds to the exponent of the local power law. Again, one can see the areas 
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ZRP forcing , wind speed U=10 m/sec 




xg/U A 2 


(a) Solid line - numerical experiment, (b) Exponent q of the mean frequency 
dashed line fit by 3.4 • (/■)§■) ° 3 . dependence on fetch x 


Figure 2 


corresponding to Kolmogorov-Zakharov index —4 and Phillips index —5. The 
value of the index to the left side from —4 plateau has the tendency to grow, 
which qualitatively corresponds to the “inverse cascade” Kolmogorov-Zakharov 
index —11/3. 

One should stop on e ~ w~ 4 asymptotics. It was observed in all our numer¬ 
ical experiments. It is Zakharov-Filonenko spectrum, which is the solution of 
equation 


Snl = 0 


(23) 


It is predicted by weak-turbulent wind-wave turbulence theory and appears 


routinely in numerical experiments 
This spectrum is confirmed by multiple ocean field 



wave tanks 


5l| and Lake George 29( measurements. 


The ’’inverse cascade” spectrum £ w ~ w 11 / 3 was also predicted by weak- 


turbulent theory 


H, HH 


fl¬ 


its 


and observed in numerical experiments 
field measurements, however, are less confident. 

In reality, nonlinear S n i term is the leading term in the ocean energy balance 
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In f 


(a) Logarithm of spectral energy den¬ 
sity as a function of logarithm of fre¬ 
quency f — - solid line. Dashed 

line - fit ~ w -4 , dash-dotted line - fit 




LnF 


(b) Local exponent of to calculated 
from the solid line Fig ldal 


Figure 3 



It consists of two parts: 


S n i = Fk - T k e k (24) 

which almost compensate each other. Otherwise, one can not explain persistent 
presence of Zakharov - Filonenko asymptotics ~ w -4 . 

Fig 0 presents relation ”(10q — 2 p) as a function of fetch x. It is in perfect 
accordance with self-similar prediction Eq. m- 

We conclude that alternative framework for HE simulation reproduces the fol¬ 
lowing analytical features of HE: 

1. Self-similar solutions with correct exponents 

2. Kolmogorov-Zakharov spectra ~ w -4 

Table [T] presents results of calculation of exponents p and q (see Eas. (fT7|) - (fT5|l ') 
for 14 different experimental observations with the last row corresponding to 
limited fetch growth numerical experiment within alternative framework. One 
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ZRP forcing , wind speed U=10 m/sec 



Figure 4: Relation (10q-2p) as a function of the fetch x. 

can see good correspondence between theoretical, experimental and numerical 
values of p and q. 

5. Tests for separation of trustworthy wind input terms from non¬ 
physical ones 

As it was already discussed, there are plenty of historically developed param- 
eterizations of wind input terms. Analysis of nonlinear properties of HE in 
the form of specific self-similar solutions and Kolmogorov-Zakharov law for di¬ 
rect energy cascade allows us to propose the set of tests, which would allow 
separation of physically justified wind-input terms S, m from non-physical ones. 
As such, we propose: 

1. Checking powers of observed energy and mean frequency dependencies 
along the fetch versus predicted by self-similar solutions. 

2. Checking the “Magic relations” Ea. (l2lTl between exponents p and q for 
observed energy and frequency dependencies along the fetch. 

3. Checking exponents of directional spectral energy dependencies versus 
Kolmogorov-Zakharov exponent —4. 

We applied such tests to the results of HE simulations which used the following 
popular wind input terms within alternative framework: 
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Experiment 

P 

q 

Babanin, Soloviev 1998 

0.89 

0.28 

Walsh et al. (1989) US coast 

1.0 

0.29 

Kahma, Calkoen (1992) unstable 

0.94 

0.28 

Kahma, Pettersson (1994) 

0.93 

0.28 

JONSWAP by Davidan (1980) 

1.0 

0.28 

JONSWAP by Phillips (1977) 

1.0 

0.25 

Kahma, Calkoen (1992) composite 

0.9 

0.27 

Kahma (1981, 1986) rapid growth 

1.0 

0.33 

Kahma (1986) average growth 

1.0 

0.33 

Donelan et al. (1992) St Claire 

1.0 

0.33 

Ross (1978), Atlantic, stable 

1.1 

0.27 

Liu, Ross (1980), Lake Michigan, unstable 

1.1 

0.27 

JONSWAP by Hasselmann et al. (1973) 

1.0 

0.33 

Mitsuyasu et al. (1971) 

1.0 

0.33 

ZRP numerics 

1.0 

0.3 


Table 1 
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(a) 


(b) 


Figure 5: Same as Figd but for Chalikov 


1 . 

2 . 

3. 

4. 


Chalikov Si n term 
Snyder Si n term 
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15 | 


301 


Hsiao-Shemdin Si n term 
WAM3 S in term [521 


m 


6. Test of Chalikov wind input term 

Fig[5] shows that total energy growth along the fetch significantly exceeds ob¬ 
served in ZRP simulation, and value of the corresponding exponent significantly 
deviates from theoretical value p = 1.0 

Dependence of mean frequency against the fetch, shown on FigEl also deviates 
from ZRP numerical results and corresponding self-similar exponent q = 0.3 
Left side of FiglTal presents directional spectrum as a function of frequency in 
logarithmic coordinates. One can see that similar to ZRP case we observe: 

1. Spectral maximum area 

2. Kolmogorov-Zakharov segment ~ uj~ 4 

3. Phillips high frequency tail ~ w -5 

Fig ITbl shows log-log derivative of the energy curve from FigjTal corresponding to 
the exponent of the local power law. Again, one can see the areas corresponding 
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(a) (b) 

Figure 6: Same as Fig[21 but for Chalikov 

to Kolmogorov-Zakharov index —4 and Phillips index —5. The value of the 
index to the left side of —4 plateau has a tendency to grow, which qualitatively 
corresponds to the Kolmogorov-Zakharov inverse cascade index —11/3. 

FigH] presents combination (10g — 2 p) as a function of fetch distance x. It is 
surprising that it is in perfect accordance with the relation Ea. (121l) . It mean that 
despite incorrect values p and q along the fetch, their combination (10g — 2 p) 
still holds in complete accordance with theoretical prediction, i.e. self-similarity 
is fulfilled locally. 

7. Test of Snyder wind input term 

Fig[9] shows that total energy growth along the fetch significantly exceeds ZRP 
case, but has the value of growth exponent close to p = 1.0 versus fetch coordi¬ 
nate x. 

Dependence of mean frequency against the fetch shown on Fig llOl is lower than 
ZRP numerical results, but has fairly close value to self-similar solution index 
q = 0.3. 

Fig lllal oresents directional spectrum as a function of frequency in logarith¬ 
mic coordinates. One can see: 

1. Spectral maximum area 
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(a) (b) 

Figure 7: Same as Fig.3, but for Chalikov S',; n 



Figure 8: ’’Magic number” lOg — 2p as a function of the fetch x for Chalikov 
wind input term. 
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(a) (b) 

Figure 9: Same as Fig[l] but for Snyder Si n 




(a) (b) 

Figure 10: Same as Figj2j but for Snyder Si n 
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(a) (b) 

Figure 11: Same as Fig 13] but for Snyder Si n 


2. Kolmogorov-Zakharov segment ~ w 4 

3. Phillips high frequency tail ~ w -5 

Fig lllbl shows log-log derivative of the energy curve from Fig lllal which cor¬ 
responds to the exponent of the local power law. Again, one can see the areas 
corresponding to Kolmogorov-Zakharov index —4 and Phillips index —5. The 
value of the index to the left side of —4 has a tendency to grow, which qualita¬ 
tively corresponds to the “inverse cascade” Kolmogorov-Zakharov index —11/3. 
Fig HU presents the combination (10g — 2 p) as the function of the fetch. Again, 
it is in perfect accordance with the theoretical relation Eq. (1711) . As in Chalikov 
case it means that despite not perfect values of p and q and wrong energy growth 
along the fetch, their combination (10g — 2 p) still holds in complete accordance 
with theoretical prediction, i.e. self-similarity is also fulfilled locally in Snyder 
case. 


8. Test of Hsiao-Shemdin wind input term 

FigHU shows that total energy growth along the fetch strongly underestimates 
ZRP simulation, and has the asymptotic value of exponent p ~ 0.5. 
Dependence of the mean frequency against the fetch shown on Fig[l4] demon¬ 
strates discrepancy with ZRP results and asymptotic value of index q ps 0.21. 
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Figure 12: Relation (lOg — 2p) as a function of the fetch x for Snyder wind input 
term. 




(a) 


(b) 


Figure 13: Same as FigQ] but for Hsiao — Shemdin Si n 
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(a) (b) 

Figure 14: Same as Figj2] but for Hsiao — Shemdin Si n 

Fig ll5l presents directional spectrum as a function of frequency in logarithmic 
coordinates. One can see: 

1. Spectral maximum area 

2. Kolmogorov-Zakharov segment ~ w -4 

3. Phillips high frequency tail ~ w -5 

Fig ll5bl shows log-log derivative of the energy curve from Fig llSal corresponding 
to the exponent of the local power law. Again, one can see the areas correspond¬ 
ing Kolmogorov-Zakharov index —4 and Phillips index —5. The value of the 
index to the left side of —4 plateau has the tendency to grow, which qualitatively 
corresponds to the “inverse cascade” Kolmogorov-Zakharov index —11/3. 

Fig HU presents combination (10g — 2p) as the function of the fetch coordinate x. 
It is in total agreement with the theoretical predictions Eq. m, which means 
that self-similarity is fulfilled locally in Hsiao — Shemdin case. 

9. Test of WAM3 wind input term 

Fig d7] shows that total energy growth along the fetch dramatically underesti¬ 
mates ZRP simulation, and has the value of exponent p asymptotically going 
to 0 versus fetch coordinate x. 
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(a) 


(b) 


Figure 15: Same as Fig[3]but for Hsio — Shemdin Si n 



Figure 16: Relation (lOg — 2 p) as a function of the fetch x for Hsiao — Shemdin 
wind input term. 
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(a) 


(b) 


Figure 17: Same as Figd but for WAM 3 Si n 


Dependence of the mean frequency against the fetch shown on FigjTB] demon¬ 
strates strong discrepancy with ZRP results and corresponding index q also 
goes to 0 asymtotically. 

Fig ll9al presents directional spectrum as a function of frequency in logarithmic 
coordinates. One can see: 

1. the spectral maximum area 

2. Kolmogorov-Zakharov segment ~ w -4 

3. Phillips high frequency tail ~ w -5 

Fig ll9bl shows log-log derivative of the energy curve from Fig ll9al corresponding 
to the exponent of the local power law. Again, one can see the areas correspond¬ 
ing Kolmogorov-Zakharov index —4 and Phillips index —5. The value of the 
index to the left side of —4 plateau has the tendency to grow, which qualitatively 
corresponds to the “inverse cascade” Kolmogorov-Zakharov index —11/3. 

Fig [20] presents combination (10g — 2 p) as the function of the fetch coordinate 
x. It is in total disagreement with the theoretical predictions. There is no any 
indication of “magic relation” Ea. (l2lT) fulfillment. 

Comparing the results, obtained for Snyder and WAM 3 wind input terms, we 
see strong discrepancies. Energy dependence on the fetch for WAM3 model 
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(a) (b) 

Figure 18: Same as Fig[2l but for WAM3 Si n 




(a) (b) 

Figure 19: Same as Fig[5]but for WAM 3 
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Figure 20: Combination (10g — 2p) as a function of the fetch x for WAM 3 wind 
input term. 

(which is esentially Snyder model with long-wave dissipation) strongly qualita¬ 
tively differs from Snyder model. When long-waves dissipation is present, the 
energy is not only much smaller, but also its dependence on the fetch stops to 
be power-one very soon, and radically differs from all other above considered 
variants of wind forcing, for which long-wave dissipation is absent. 

10. Conclusion 

We are offering alternative framework for numerical simulation of HE. Being 
supplied with ZRP wind input term, such approach reproduces the results of 
more than a dozen of experimental observations. 

We also performed numerical simulations of HE for four another historically 
well-known wind input terms within the same alternative framework. They 
demonstrated the results deviating from ZRP simulation. 

To classify the results of the above simulations we applied the set of nonlinear 
tests to different kinds of wind input terms, and here is the conclusion: 

1. ZRP forcing term perfectly corresponds to theoretically predicted results 
like Kolmogorov-Zakharov spectrum ~ , self-similar solutions for 
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Experiment 

p-test 

q -test 

A'Z-spectrum 

Magic relation 

Energy growth 

ZRP 

YES 

YES 

YES 

YES 

YES 

Chalikov 

NO 

NO 

YES 

YES 

NO 

Snyder 

S3 

S3 

YES 

YES 

NO 

Hsiao — Shemdin 

NO 

NO 

YES 

YES 

NO 

WAM3 

NO 

NO 

YES 

NO 

NO 


Table 2 

energy and frequency with exponents p = 1 and q = 0.3 correspond¬ 
ingly, “Magic relation” lOp — 2q = 1 and reproduces more than a dozen of 
field experiments. Therefore, it can serve as the benchmark. 

2. All wind input terms pass the test for presence of Kolnrogorov-Zakharov 
law e ~ uj~ 4 . This means that effects of nonlinearity are so strong, that 
presumably no variation of the wind input term parameterization can 
suppress it. 

3. Chalikov and Hsiao — Shemdin cases fail p — and q— tests, but pass 
“Magic relation” (quasi-self-similarity) test. 

4. Snyder case “approximately” passes p—, q— and “Magic relation” tests. 

5. W AM 3 case fails to pass all exept KZ— spectrum test. 

6. None of the wind-input parameterization, except ZRP one, can correctly 
reproduce experimentally observed limited fetch growth. 

In summary, the nonlinearity influence is so robust in the dynamics of HE that 
one can’t “spoil” Kolnrogorov-Zakharov law ~ cv~ 4 for any tested wind input 
term S\ n . Self-similarity tests like p — and q— tests are the most sophisticated 
between suggested ones. And the “magic relation” test is probably somewhere 
in-between versus detection of the “quality” of particular wind input term. 

The summary of the tests is presented in Table [2] 
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